ANL-HEP-PR-03-005-Rev 
EFI-03-04-Rev 
hep-ph/0302094 
March 2003 



NEW PHYSICS CONTRIBUTIONS TO THE B -> (f)Ks DECAY 1 

OH ! Cheng- Wei Chiang 

High Energy Physics Division, Argonne National Laboratory, Argonne, IL 6044 
1 and 
?h . Enrico Fermi Institute, University of Chicago, Chicago, IL 60637 

<^ Jonathan L. Rosner 

Enrico Fermi Institute and Department of Physics 
University of Chicago, Chicago, IL 60637 

m 
> 

On 

Recent measurements of the time- dependent CP asymmetry of the B — > 4>Ks 
decay give results whose central values differ from standard model expecta- 
tions. It is shown how such data can be used to identify new physics contribu- 
tions in a model-independent manner. In general, a sizeable new amplitude 
JpH. with nontrivial weak and strong phases would be required to explain current 

Qh! data. Improvement in the quality of data will allow one to form a more 

J|h ' definite conclusion. 



X 



ABSTRACT 



Submitted to Phys. Rev. D. 



1 



I. INTRODUCTION 



It has been suggested to look for discrepancies among the time-dependent CP asym- 
metries of different B decay modes as a means to detect new physics [H 121 El IH El E] • 
Since the B — > J/ipKs decay is a tree-dominated process in the standard model (SM), 
its CP asymmetry Sj/^k s is believed to be less affected by new physics and to give infor- 
mation on sin 2(3. Although the CP asymmetry of the B —>■ <pKs mode is also expected 
to give the same sin 2(3 within the SM, this process is, however, particularly sensitive to 
new physics contributions because it is a purely penguin loop-mediated process in the 
SM. The SM pollution from a small w-penguin with the weak phase 7 has been studied 
in Ref. [7] and it is found that the deviation of S<j,K s from sin 2(3 is of 0{\ 2 ) ~ 5%, where 
A ~ (9(0.2) is a parameter close in magnitude to the Wolfenstein parameter A ~ 0.22 8J. 
Therefore, a large deviation of S^Ks from its SM prediction would signal contributions 
from physics beyond the SM. 

As argued by Fleischer and Mannel [I] , even if one ignores rescattering effects, contri- 
butions from new physics at a TeV scale to AI = operators could be of the same order 
as the SM ones, while new AI = 1 operators are suppressed by A. With rescattering 
effects taken into account, both the new AI = 1 operators and the SM pollution will 
be enhanced by about A. In any case, both of the AI = 0, 1 operators from the new 
TeV-scale physics can be more significant in comparison with the above-mentioned SM 
pollution. 

The world average of sin 2(3 as measured from the golden mode B — > J/ipKs, 
sin 2(3 = 0.734 ± 0.054 jH], agrees well with constraints obtained from other experi- 
ments. Recently both the BaBar and Belle groups have also reported measurements 
of time-dependent CP asymmetries in the B — > (pK s decay. S^Ks (the coefficient of 
sin Ami in flavor-tagged decays) is found to be about 2.7a away from Sj/,pK s , while 
•A-4>k s (the coefficient of cos Ami) is la away from 0. If this situation continues as the 
data precision improves, it would be of interest to know the magnitude and phase of 
possible new physics contributions to the <pKs mode. 

Instead of separating the new physics contributions into AI = and AI = 1 parts as 
done in Ref. j3], we will simplify the discussion by considering the combined amplitude 
from such effects along with the smaller SM pollution amplitude. This enables us to 
obtain useful information from the three observables <S</>k s , ^k s , and the ratio R be- 
tween the sum of squared amplitudes extracted from the measured B°(B ) — > <f)K°(K ) 
branching ratios and a "standard" squared amplitude, such as the SM predicted value or 
experimentally measured B ± — > K*^ branching ratio. The algebraic structure of the 
problem then becomes very similar to that studied by several authors JTU] for B — ► tttt. 
We try to find in a model-independent way the allowed magnitude and phases of the 
new amplitude and some generic properties associated with it. Such an analysis is useful 
in helping us narrow down new physics models jTTj consistent with observed data. 

The paper is organized as follows. Sec. II introduces a decomposition of decay am- 
plitudes in terms of topological contributions. The formalism for time-dependent CP 
asymmetries is discussed in Sec. III. We present numerical analyses for two separate 
cases of new physics in Sees. IV and V. In Sec. VI, we summarize our results. 
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II. TOPOLOGICAL AMPLITUDE ANALYSIS 



In the framework of the SM, both the <J)K° and <j)K modes receive important con- 
tributions from QCD and EW penguin graphs, with the former having a dominant 
effect. Useful information about the QCD penguin contribution can be obtained from 
the K*°7T decay mode using fiavor-5'C/(3) symmetry jT^j- It should be noted that a tiny 
annihilation diagram also exists in both the (pK^ and K* ^ decay modes. From the 
arguments of both dynamical suppression and the fact that no asymmetry is observed 
between the K* ^ modes, we shall ignore the annihilation amplitude in these charged 
decays. In this case, both the neutral and charged <pK modes have the same decay 
amplitudes. We will then average over the branching ratios of these two sets of modes 
using their associated errors as the weights for our analysis. 

Let's write down the amplitudes of the relevant modes in terms of independent 
topological components as [T31 Ej 

A(K*\ + ) = pe l ^ SM+&p) , (2) 

in the SM. In the above two equations, the p part denotes the QCD penguin contribution 
which also contains a negligible color-suppressed EW penguin amplitude, and the s part 
denotes the EW penguin contribution along with a small flavor-S'?7(3)-singlet amplitude, 
as expected from the OZI rule. The variables p and s are absolute values of the respective 
amplitudes and therefore are non-negative by definition. The weak phase <psu satisfying 
e -2i<t>sM — V t bV t * s / (y t * b Vt s ) is the same for both the p and s parts. Finally, 5 P and 5 S are 
the associated strong phases. Note that to simplify the notation given in Ref. |12j . we 
omit from these amplitudes the subscript P indicating that the spectator quark ends 
up in the pseudoscalar meson in the final state and the prime that denotes AS = 1 
transitions. 

It should be noted that we explicitly assume flavor SU(3) symmetry in Eqs. (JTJ) and 
(J21) in order to relate the amplitude for the K*°tt + mode to the penguin part in the <f)K 
process in later analysis. The SU(3)f breaking effect will be characterized by the factor 
[f$F B ^ K (m 2 ^)]/[f K *F B ^(m 2 K ,)} ~ 1.2. We will simply treat this extra factor as 1 in 
our analysis. 

New physics can give rise to new operators that contribute to the decays of the 
above processes. We will distinguish two cases in later discussions: (i) only the <pK 
modes receive the new contributions while the i^*°7r ± modes are purely SM processes; 
and (ii) both types of decay modes receive the same contributions from new physics. 
Case (i) could happen, for example, when new physics enters the h — > s EW penguin 
only. In this case, we add an extra amplitude 77, e l( ^ n+<5n ) to Eq. (0). Case (ii) could 
happen when new physics modifies the b —>■ s QCD penguin. In that case we add the 
new amplitude to both Eq. and Eq. (j2J). In general, the new amplitude 77, e l( ^ n+<5n ) is 
a combination of AI = and AI = 1 ones that may contribute at different strengths [3] . 
A more careful job can in principle be done by separating the new amplitude into those 
with different isospins and studying new physics contributions in each piece. However, 
one would find that there are not enough observables among the decay modes to solve 
for all the parameters in the amplitudes. 



3 



III. TIME-DEPENDENT CP ASYMMETRIES 

In this section, we review the general analysis of time- dependent CP asymmetry of 
pure B° and 5° decays into a CP eigenstate fcp- Let's define the asymmetry as 



, t) = F ffiU*) - fcp) ~ r (gUg ^ fcp) 

afcp r « ys (t) - /cp) + r (B° hys (t) - / CP ) ' 



In our case, fcp — 4>Ks- Denote 



where ?70k s = — 1 is the CP eigenvalue of the (f>Ks state, 



are factors that account for the mixing effects in neutral B and K meson systems, 
respectively, and 

A((j)K ) = (K°\H\B°) =ae iW « +y + be l( ^ +Sb) , (6) 

A((fK°) = (K°\n\B°) =ae i( -* a+6a) + be l (- <t > b+Sb) , (7) 

where a, b are chosen to be positive, 4>a,b G {— 7r, iy} and 5 aj6 G {0, 2ir} are the associated 
weak and strong phases, respectively. The above amplitudes are invariant under the 
transformations <p at b — > a) ft±m7r, 5 a ^ — > 5 aj 6=Fm7r and ai 5 — > ai &±m7r, 5 a 6 — > 5 aib ±m7r 
for m G Z. Here the separation of the total amplitude into two parts is done in accord 
with the nature of the problem. The ratio of the amplitudes in Eqs. (6) and (7) is then 



A((j)K ) 1 + re-^+ & ) ' { ' 

where 

r = b/a>0 , = a - 05 , and <5 = 5 a - 5 6 . (9) 
The CP asymmetry can then be written as 

a^Ksit) = AcpKg cos(AMt) + S<p Ks sin(AMt) , (10) 

where AM is the mass difference between the two physical B meson states, and 

A<t>K s = l/ yg L - -i , (11) 
2 ImA^ 
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New physics will affect the CP asymmetry observables through the parameter X^k s - 
Therefore, it may come in at two places: the mixing matrix and/or the decay amplitudes. 
As emphasized in Ref. pQ, new physics effects on the mixing part will be universal and 
do not change the SM predicted pattern of CP asymmetries in different modes; their 
effects on the decay amplitudes, however, are non- universal so that the CP asymmetries 
can vary from channel to channel. Since current sin 2(3 measurements from other decay 
modes, such as J/ipKs, rj'K s , etc, seem to agree with one another and with the unitarity 
triangle constraints obtained from other processes pretty well, it is plausible to assume 
that any strange behavior in the (fiKg mode is mostly due to new physics contributions 
in the amplitudes. In this case, we will use the SM mixing factors in Eq. (j3J) and obtain 

1 _i_ rp i((f>-6) 

^ = - e " m - 1+re -,(^) . < 13 > 

where 

e-MA* = Q Q e -^ a (14) 

Within the SM, <fi a ~ tt and one obtains the effective weak phase [3 e g coinciding with f3 
in the unitarity triangle. However, if new physics modifies the phase (f) a , then /3 e g will 
in general differ from what the SM expects. If one writes 4> a = 4>sm + 4> with (fisM being 
the phase expected in the SM and (fi being the deviation, then j3 e s = cfisu + (fi (mod tt). 

We will be exclusively dealing with the decays of a B meson into a final state with 
one pesudoscalar meson (P) and one vector meson (V). The invariant amplitude A of 
such a process is conventionally related to its partial width in the following way: 

T(B -> PV) = ^Xt\A{B -> PV)\\ (15) 



where p* is the 3-momentum of each final particle in the rest frame of the B meson, and 
is the mass of the decaying B meson. Note that p* is raised to its third power to 
appropriately account for the P-wave kinematic factor. 

IV. NEW PHYSICS ONLY IN (fiK SYSTEM 

In this section, we will discuss the case when new physics only enters the (fiK system 
but not the K*tt system. Since the p and s parts of the (fiK decay amplitudes have the 
same weak phase (fisM, we can combine them into a single v part and write, including 
the new physics part, 



A((pK u ) = v e ^ sM > + n e' 1 ^ = v e l[cpSM+dv) 1 + r e ~ Wj , (16) 

where r = n/v, <fi = 4>sm — (fin, and 5 = 5 V — 5 n . Since here we assume that the K* ^ 
modes are not affected by the new physics, we can use them to obtain reliable information 
on p, the magnitude of the QCD penguin. The effective Hamiltonian approach indicates 
that there is a small relative strong phase between the p and s amplitudes in the SM 
[T51 116j. The relative strong phase obtained in this approach comes purely from short- 
distance physics [T7]. In general, there are nonperturbative strong phases from soft 
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gluon exchanges in the final-state particles that may be different between the two types of 
penguin diagrams. For simplicity and definiteness, we will take 8 V = 5 P ~ 8 S — 7r ~ since 
the overall strong phase will not matter, and consider maximal destructive interference 
between the QCD and EW penguins, in accord with the effective Hamiltonian analysis. 
The additional — 7r in the above strong phase relation does not really have a strong 
interaction origin but simply comes from the charge coupling of the final-state s quark 
with the Z boson in the EW penguin. Therefore, the s part has a 180° phase from 
the p part within the SM. Under this assumption, v = p — s can be computed once we 
know the prediction of the ratio s/p in the SM. We will also mention consequences of 
imperfect destructive interference between these two types of amplitudes. 

A. Observables 

One observable in the 0A system is the ratio 

The numerator in the above definition is the sum of measured branching ratios of B° — > 
4>K° and B — > <f)K . As described before, we will actually take the weighted average of 
the neutral and charged modes. The denominator is the theoretical prediction for the 
same branching ratio sum within the SM. In terms of R and Eq. (|13jh we obtain 

RS+Kg = sin2/? + 2rcos5sin(2/3 - 0) +r 2 sin2(/3 - 0) , (18) 
RA^Ks = 2 r sin sin 5 . (19) 

Now we have three observables R, S^k S i and A$k s that allow us to solve for the three 
parameters r, 0, and 5. As the value of R may vary owing to the interference between 
p and s, we will estimate the SM contribution and also search the allowed parameter 
space by varying R over a reasonable range. 

The self-tagging modes 5 ± — > (pK^ can provide additional statistical power to the 
determination of A^Ks if we assume that A(B + — > (f>K + ) = A(B° — ► 0A°) as we have 
done above. In that case one finds just 

for the time-integrated CP rate asymmetry. The BaBar Collaboration JH] has recently 
reported Aqp = 0.039 ± 0.086 ± 0.011. We shall not use this value in our averages but 
in principle it can greatly reduce the error on A</>k 3 - 

B. Numerical studies 

In this subsection, we will use the measured values of A^Ks, S<i>k s , and R with some 
theoretical input from the SM to find the allowed ranges of r, 0, and 5. Solving for r in 
Eq. ()17|) in terms of R, 0, and 5, one obtains two solutions 



r 1 = — cos cos 5 — J cos 2 cos 2 5 + R — 1 , (Solution I) 



r 2 = — cos cos 5 + a/ cos 2 cos 2 5 + R — 1 . (Solution II) (21) 
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Table I: Experimental input of measured branching fractions. 



(xlO- 6 ) 








B(<f>K+) 




CLEO 


5.4±|;^ ± 0.7 (< 12.3) [E] 


5.51^ ±0.6 [IH] 


7.6ii:g ± i.6 [U] 


BaBar 


7.6tJ-j ± 0.5 DSI 


io.o±8s ± o.5 [IH! 


15.5 ±3.4 ±1.8 [21] 


Belle 


in n+ L9 + - 9 

1U - U -1. 7-1.3 


22. 




10.7±1.0i?| |22] 


1 q /( +4.2+2.1+3.5 ryfri 
iy-4_3. 9 _2.i-6.8 


Average 


7.98 ±1.07 




9.51 ±0.78 


12.3 ±2.5 



Table II: Experimental input of measured CP asymmetries for the B — > 4>Ks mode. 



Quantity 


BaBar 


24 


Belle 25J 


Average 


S 


-0.18 ±0.51 ±0.07 


-0.73 ±0.64 ±0.22 


-0.38 ±0.41 


A 


0.80 ±0.38 


±0.12 


-0.56 ±0.41 ±0.16 


0.19 ±0.30 



First, it is seen that Solution I is not allowed for R > 1 because t*i has to be positive. 
Therefore, we see that r\ < 1. 

In the effective Hamiltonian approach, the EW penguin is found to be of consider- 
able importance, with the ratio \s/p\ predicted to be between 10% and 11% using the 
results given in Ref. ^H]- On the other hand, the B + — > K*°tt + decay mode involves 
only p, ignoring a small annihilation diagram that also contributes to the <pK + mode. 
Using its branching ratio, we obtain \p\ = (1.42 ± 0.14) x 10 -8 . Combining the above 
results and assuming maximal destructive interference between p and s, the SM predicts 
\A SM ((j)K)\ = (1.27 ±0.13) x 10~ 8 . To improve the statistics, we take the weighted 
average for the branching ratios of the neutral and charged <j)K modes as given in Table 
[Hand obtain \A cxp (<pK)\ = (1.27 ± 0.04) x 10 -8 after removing the kinematic factors. 
Therefore, we obtain an estimate of 



R 



A exp ((pK) 



A SM {<pK) 



0.99 ±0.21. (22) 



If a nontrivial relative strong phase exists between p and s, the central value of the 
resulting R will become smaller. In the case of maximal constructive interference between 
p and s, R could be as low as 0.5. 

We will use the CP asymmetry observables measured by the BaBar and Belle groups 
[2*4^ 123] as given in Table ITT1 for our analysis. Replacing r in Eqs. ()18|) and (fTT?[) by one of 
the above solutions, it is then possible to find on the <j)-5 plane regions that are consistent 
with the measured values A^ Ks = 0.19 ± 0.30, S (j)Ks = -0.38 ± 0.41 [21 [23, and the 
additional requirement that r > by definition. The fact that A^,k s is negative at 
the la level gives the following possibilities: (i) — tt < < and < 8 < tt; and (ii) 
< < 7T and tt < 5 < 2tt. 

In Fig. ^ we only show a set of representative solutions in the range — tt < <fi < 0, 
< 5 < it for R = 0.8, 1.0 and 1.2. It is noticed that Solution I does not exist when 
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71 



CO 71/2 




CO 71/2 




CO Tl/2 




-71 -71/2 71/2 



(c) 

Figure 1: The allowed regions in the <p-5 plane for Solution II with (a) R = 0.8, (b) R = 
1.0, and (c) R = 1.2. Here we only show the allowed regions in the range — ir < <fi < it, 
< 8 < n. Other regions can be obtained by <fi — > <fi ± n and 5 — > 5 ± n. Solution I 
is not allowed for R > 0.8 and, therefore, no corresponding plots are shown here. For 
R = 1 S^Ks an d A<f,K s are unchanged under <5 — ^ 7r — <5 and — > <fi + | [Fig. 1(b)]. 

8 




R > 0.8. Therefore, we only show those for Solution II. As shown in Section III, solutions 
in other regions on the <p-5 plane can be obtained by the translations <p —>■ <p ± n and 
5 — > 5 ± 7r. As R increases, the allowed regions of Solution II become larger. 

In Fig. El we plot the allowed ranges of r for < R < 1.4. The dark region in plot 
(a) corresponds to Solution I, and that in plot (b) to Solution II. It is seen from the 
plots that to satisfy the constraints of measured data, r has to be at least about 0.4 for 
either solution. This corresponds to a new physics amplitude with a magnitude of at 
least about 0.45 x 10 -8 . It is also found that for Solution I, R has to be less than about 
0.8. Therefore, the current value of R favors Solution II. If we take the value R — 1, 
Solution II has a wide range for r: 0.40 ^ r < 0.90 and 1.05 ^ r < 1.96. 

The fact that r has to be greater than a minimum can be readily understood. Should 
r be too small, then new physics [the n part in Eq. (fTU|) ] does not have enough weight to 
change RS<f>K s from that extracted from the J/ipKs mode to the measured one as the 
modification is of 0(r) according to Eq. (|18|). As mentioned in the beginning, r ~ 0(1) 
means that the new amplitude has the same order of size as the SM contribution. This 
would point to the possibility of new physics at the TeV scale or below. 

The sensitivity of S^Ks an d A^,k s to the weak and strong phases and 5 for a value of 
R close to the central one is illustrated in Fig. El Here each curve for a given </> intersects 
the axis A$k s = at either 5 = or 5 = it, while curves with A$k s < are related 
to those with A$k s > and the same value of S^Ks by the transformation (f> — > n + (ft, 
5 —>■ 7i — 5. The plotted cross shows the present status of the data summarized in Table 
II. The ranges of 5 and <fi are restricted in general by the requirement that the argument 
cos 2 (f> cos 2 S + R — 1 of the square roots in Eqs. (J2Tj) be non-negative. Because of the 
special value of R, we are able to draw curves for essentially any value of by varying 
5. Therefore, the constraints come merely from S^Ks and A^k s - 
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(3 = 23.6°, R = 1.0, large r 




-1.0 -0.5 0.0 0.5 1.0 



S 

Figure 3: The curves traced in the S^Ks'^Ks plane by varying the relative strong phase 
S between and ir for fixed values of 0. The plot is for (5 = 23.6° and R = 1 with r 
chosen according to Solution II; no solution I exists for R = 1. Curves are labeled by 
values of (dashed: < 0; solid: > 0) in degrees. Squares and diamonds correspond 
to values of 5 = or n. The point at S^Ks = 0.734, A^k s = corresponds to = 0, ±7r 
for all 5. The plotted data point is the average quoted in Table II. 
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TC 



tO Till 







-71 



-7C/2 







7C/2 



Figure 4: The allowed regions on the <p-5 plane for R — 1.0 using Solutions II and with 
a factor of 3 improvement in A^Ks an d S<t>K s - 



Assuming the central values of S^Ks and A$k s stay the same in future experiments 
but the errors are improved by a factor of 3, we find that the allowed regions become 
smaller. This is demonstrated in Fig. 0J 

We find that the variation of sin 2/3 within its experimental range makes little dif- 
ference in the allowed solutions. The general behavior and regions presented in Figs. H 
and El remain the same. 

V. NEW PHYSICS IN (j)K and K*n SYSTEMS 

In this section, we consider the situation of new physics entering both the <f>K and 
K*n systems. Since the new physics contribution ne % ^ n+&n ^ is to be added to both 
Eqs. (HJ) and (j2J), we combine the p and n parts into a single q part as follows 



A(K*\ + ) 

i(4>q+S q ) 



qe 



j+S q ) 



+ se 



i(4>SM+Ss) 



qe 



l((f>q+S q ) 



1+r'e 



-i(<t>+6) 



qe i(4> 9 +S q ) ; 



"> SM , and 5 



(23) 
(24) 

S n — S s . 



q YSM i " — "q 

~ 7r, one thus should use 



where qe i ^ q+Sq ^ = p e^sM+Sp) _|_ n e i{<j> n +& n ) ^ r > _ s/q, 
As mentioned in the previous section, <J)sm — 71 and 5, 
(</> g , S q ) ~ (0 + 7r, 5 + 7r) to obtain the weak and strong phases associated with the q part 
when interpreting our following plots drawn on the (j)-5 plane. 



A. Observables 



In this case, we use the observable 



-ftO 



R! 



\A cx p(B+ -> A^*°7T+)| 2 + \A cx p(B- 
1 + 2r' cos (h cos 5 + r' 2 , 



0Z°)| 2 



K 7T 



(25) 
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Note that in spite of the similarity in the forms between R' and R defined in the previous 
section, they are actually very different. Using Eq. ()25|). we have 



R'S^Kg = sm2p eS + 2r'cosSsm(2l3 e{i -(f))+r' 2 sm2(P e{i -(f)) 

= sin2(/3 + 0) + 2r'cos5sin(2/3 + 0) + r' 2 sin 2/3 , (26) 
R' A<f>K s = 2 r' sin sin 5 , (27) 

where /3 e fi = ft + is used. Here one quickly realizes that we also have only three 
parameters, r', 0, and S for which to solve. 

As in the previous case, the self-tagging rate asymmetry for — > 0fC ± provides 
additional statistical power for the measurement of A^,k s , since 

A CP = 2r' sin sin . (28) 



B. Numerical studies 

Solving r' in Eq. (|25|) in terms of R f , 0, and 5, one obtains two solutions 



r[ = — cos0 cos 5 — ycos 2 cos 2 5 + R' — 1 , (Solution I) 

= — cos cos 5 + J cos 2 cos 2 5 + i?' — 1 . (Solution II) (29) 

First, as in the previous case, Solution I is not allowed for R' > 1 because r[ has to be 
positive. Therefore, we see that r[ < 1. 

We extract the amplitude for the K*tt mode from Table to be \A exp (K*7r) \ = 
(1.42±0.14) x 1CT 8 . Combined with the weighted average of the <j)K mode amplitude size 
given in the previous section, |„4 exp (0iT)| = (1.25±0.05) x 1(T 8 , we find R' = 0.79±0.17. 
In Fig. El we show a set of representative solutions in the range —it <0<O, 0<5<7t. 
We take R' = 0.6,0.8 and 1.0. It is seen that both solutions have two allowed regions 
except for Solution I at R' = 1.0. As R' increases, the allowed regions become larger for 
Solution II. 

In Fig. we also draw the allowed ranges of r' for < R' < 1.2. The dark region 
in plot (a) corresponds to Solution I, and that in plot (b) to Solution II. It is seen from 
the plots that to satisfy the constraints of measured data, r' can go down to almost 
for R' ~ 1. If we take the value R' = 0.8, Solution I has 0.16 < r' < 0.44 while Solution 
II has a wider range 0.47 < r' < 0.97 and 1.20 < r' < 1.76. In the standard model, one 
expects r' ^ 0.1, = ±7r, 5 ~ ±7r in accord with the expected contribution (mentioned 
previously) of the electroweak penguin amplitude. 

If it turns out that r'< 1, that means the n part and the p part interfere to give 
an amplitude larger in size, agreeing with the fact that R' will be about 1. However, if 
r' ~ 0(1), then there is a cancellation between n and p such that the combined amplitude 
of the two becomes comparable to the S"?7(3)i?-singlet amplitude. Either situation would 
tell us the new physics contribution is important. 

The sensitivity of S^Ks an d A^k s t° the weak and strong phases and S in Solutions 
I and II for the case in which new physics enters through the penguin amplitude into 
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lO 71/2 



I » 



-71/2 7t/2 

4> 



(a) 



tO 71/2 



-71/2 Tt/2 

4> 



(c) 



lO 71/2 



lO 71/2 



-71 -71/2 71/2 

4> 



-71/2 



71/ 2 



(b) 



(d) 



-71 -71/2 71/2 



(e) 

Figure 5: The allowed regions on the <p-5 plane for R' = 0.6, 0.8 and 1.0, using Solutions 
I (plots (a) and (b) in the left column) and II (plots (c), (d), and (e) in the right column). 
Here we only show the allowed regions in the range —it < < 0, < 5 < tt. The other 
regions can be obtained by shifting: <fi — ■> <fi + 7r and 5 — > 5 + it. Solution I is not allowed 
for Rl — 1. In this case Fig. 5(e) displays an additional symmetry of the solution under 

0-^±7T-0. 
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2.5 



1.5 



0.2 0.4 0.6 0.8 1 1.2 
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Figure 6: The allowed range of r' for specific values of R', using (a) Solution I and (b) 
Solution II. 



both B — > cf)K and B + — > K*°ir + is illustrated in Fig. [7J As in Fig. [21 each curve for 
a given intersects the axis A.$k s — at either 5 = or 5 = it, while curves with 
•A^Kg < are related to those with A$k s > and the same value of Sct>K s by — > tt + 0, 
5 — > 7r — 5. Contrary to Fig. El not all values of and 5 are allowed for making a curve 
with a given value of R' < 1. 

It is interesting to notice that for the large-r' solution (Solution II), the curves for 
values of <p and ±7r — (f> overlap, leading to the appearance of continuity. This is because 
both curves are part of a common ellipse, obtained by solving Eqs. (J2S1), {23), and (fTTf) : 

=cos 2 + J R'-l , (30) 

where r x = 2 sin cos 2/5, and r y = 2sin0. These ellipses have their centers at coor- 
dinates ([cos 20 + R' — 1] sm2f3/R', 0). One immediately sees that the above elliptic 
equation is invariant under the transformation — >• n — <p. Since no explicit choice of 
solutions of r in Eq. (}2T?|) is made for deriving Eq. (j3*Uj) . it is valid for either solution. 
This is why each curve associated with in Fig. [7(a) is actually a portion of the cor- 
responding curve associated with n — in Fig. [7(b). The curves for both solutions are 
truncated (although not seen in the plot for Solution II because of the overlap) because 
of the conditions cos 2 cos 2 5 + R' — 1 > 0, r[ > 0. 

Again, we find that if the experimental precision of S^k s and A$k s can be improved 
by a factor of 3 with their current central values, then the allowed regions are consid- 
erably restricted. The variation of sin 2(3 within its experimental range also does not 
affect the general behavior and regions presented in this section. 

VI. SUMMARY 



R'S 4>Ks - (cos 20 + R' - 1) sin 2(5 



+ 
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(3 = 23.6°, R' = 0.8, small r' (3 = 23.6°, R' = 0.8, large r' 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



s s 

(a) (b) 

Figure 7: Same as Fig. 3 but for new physics entering through the penguin penguin 
amplitude into both B — > <pK and B + — > K* tt + . The left plot (a) is for Solution I and 
the right plot (b) for Solution II. In both plots, (3 = 23.6° and R' = 0.8. 



We have shown how to estimate the magnitude and weak and strong phases of any 
new physics contribution which might account for the deviation of the CP asymmetry 
parameters in B — > <f)K° from their standard-model values. We find that it is useful 
to compare the overall rate for this process either with that predicted in the standard 
model (a ratio R) or with that for B + — > K*°tt + (a ratio R'), which is expected to be 
dominated by the penguin amplitude in the standard model. 

It is observed from our analysis that an amplitude with considerable size, nontrivial 
weak and/or strong phases from new physics is required to fit the current experimental 
results. For example, in the case that new physics contributes to the <pK modes but not 
the K*n mode, the ratio of the new amplitude to the SM contribution has to be ^ 0.4, 
independent of the value of R. 

Current experimental data indicate that the size r or r' of a new physics amplitude 
relative to that of the standard model could well be of 0(1) for a wide range of R or R'. 
Such 0(1) parameters could indicate new physics at about the TeV scale. Of course, 
these extra contributions would need the right strong and weak phases in order to explain 
the current data. Considerable refinement of the rate and asymmetry measurements in 
B — > (pK° is necessary before the amplitude can be pinpointed satisfactorily, however. 
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